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Abstract. We present first results on event-by-event elliptic flow fluctuations in 
Au+Au collisions at ^/snn = 200 GeV obtained with the PHOBOS detector. Over 
the measured range in centrality, large relative fluctuations of 40-50% are found. The 
elliptic flow fluctuations are well described as being proportional to fluctuations in the 
shape of the initial collision region, as estimated event-by-event with the participant 
eccentricity using Glauber Monte Carlo. 
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1. Introduction 

Elliptic flow (t> 2 ) is one of the key observables in the understanding of the dynamics 
of heavy ion collisions. In the collision of two nuclei with finite impact parameter, the 
overlap region in the plane transverse to the beam direction is spatially asymmetric. 
Observation of an azimuthal anisotropy in the momentum distributions of the detected 
particles relative to the reaction plane is direct evidence of interactions between the 
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initially produced particles [1]. These interactions must occur at relatively early times, 
since the expansion of the source — even if uniform — would wash out the effect of the 
initial asymmetry. Quantitatively, the connection between the initial and final-state 
anisotropy is exploited by hydrodynamical models that relate a given initial source 
shape to the distribution of produced particles. In such calculations, one typically 
uses smooth, event-averaged, initial conditions. However, event-by-event fluctuations 
in the shape of the initial interaction region must not be neglected. As a means to 
quantify the effect of initial-state eccentricity fluctuations, we have introduced the 
"participant eccentricity", e part . It unifies the eccentricity-scaled elliptic flow across 
Cu+Cu and Au+Au collisions [2, 3]. It accounts for nucleon-position fluctuations in the 
participating nucleon distributions by calculating the eccentricity, event-by-event, with 
respect to the principal axes of the overlap ellipse in a Monte Carlo Glauber (MCG) 
simulation. In a hydrodynamical scenario, such fluctuations in the shape of the initial 
collision region would lead naturally to corresponding fluctuations in the elliptic flow 
signal. Note that the measurement of v 2 fluctuations as a probe of early stage dynamics 
of heavy ion collisions has been suggested earlier [4] , although without considering the 
possibility of initial-state eccentricity fluctuations. In order to experimentally determine 
elliptic flow fluctuations, we have developed an analysis technique for the extraction of 
dynamical flow fluctuations [5]. Our method, summarized below, relies crucially on the 
large pseudo-rapidity coverage of the PHOBOS multiplicity array, which detects a very 
large fraction of the produced particles [6]. The results as a function of centrality in 
v /5nn = 200 GeV Au+Au collisions at mid-rapidity are for the first time presented in 
these conference proceedings. 

2. Elliptic flow fluctuations analysis technique and results 

In order to disentangle known (mostly statistical) from unknown (dynamical) contri- 
butions to the measured flow fluctuations, a detailed knowledge of the detector response 
is required. The response function, K{vf >s 1 v 2l n) 1 is defined as the distribution of the 
event-by-event observed elliptic flow, t>2 bs , for events with constant input flow value, v 2 , 
and multiplicity, n. It is designed to account for detector deficiencies, as well as for 
multiplicity and finite-number fluctuations. Let f{v 2 ) be the true v 2 distribution for a 
set of events in a given centrality class, and g{v 2 hs ) the corresponding distribution of 
the observed v 2 , i>2 bs - The true and observed distributions are related by 



where N(n) is the multiplicity distribution of the events in the given set of events. Our 
analysis method consists of three steps: 

i) Measurement of the v 2 hs distribution, g(v 2 hs ): The event-by-event measurement 
utilizes a two-parameter maximum likelihood fit to the obtained hit information in 
an event. The probability distribution function (PDF) of a particle to be emitted 
in the direction (//,</>) for an event with mid-rapidity flow amplitude v 2 = v 2 (0) and 




(1) 
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reaction plane angle 0o is defined as P(n,(f)]V2,(po) = s 1 [1 + 2t> 2 (r/) cos(20 — 20 o )]- 
The pseudorapidity dependence, v 2 (r]), is parametrized with a triangular shape, 

v^iv) = v 2 (1 — or alternatively with a trapezoidal shape, v t 2 ra ' v (7]) = {a^^eise ' 
Both parametrizations yield good descriptions of the previously measured (mean) U2O7) 
shapes [7]. The normalization (s _1 ) is constructed such that the PDF folded with the 
PHOBOS acceptance (A cc ) yields the same value for different sets of parameters (v 2 , (f>o), 
s(v2, 0o; v) = fv-Arj/2 Acdjf ■> 0') [1+2^2(7/) cos(20'— 20 o )] d0' dn' , in practice constructed 
in small bins of An. The desired value is given by the pair (v 2 bs , 4>Q hs ) that is "most 
likely" for the given event, found by maximizing nr=i s ^iVii <fe u 2> 0o) as a function of 
t>2 and 0o for all hits in the event. 

ii) Construction of the kernel, K(v% hs , t>2, n): The response function K(v% hs , t>2, n) is 
determined by running the event-by-event analysis, as described above, on modified 
HIJING events with known values of t>2 and n. Flow of constant magnitude (ug) 
with a flat reaction plane distribution (0o) is introduced into the generated HIJING 
events. This is achieved by redistributing the resulting particles in each event in <fr 
according to 1 + 2v 2 {vj) cos(20 — 20 o ), using their generated r\ positions. The modified 
HIJING events are fed through GEANT to obtain the PHOBOS detector response. For 
each parametrization of V2{i]), triangular and trapezoidal, used in the event-by-event 
measurement, a corresponding kernel, K tri and K tra,p , is constructed. To obtain the 
kernel in bins of t>2 and n with enough precision would require on the order of 100 million 
MC events. Instead the kernel can be parametrized, allowing us to use about 1.5% of 
that statistics to reach the required precision. The response of a perfect detector can 
be determined as a function of event multiplicity as described in Ref. [8] (see Eq. (A13) 
with a — > t>2 bs , & —>■ V2 and M — > n). In practice, however, it turns out that t>2 is 
suppressed, with the suppression dependent on n, and that the resolution (cr) has a 
constant background contribution. With = (An + B) t>2 and a = C/^/n + D, this 
leads to 

= t exp ( J<)\\^ ) h (J^L\ (2) 



a 2 \ 2a 2 I \ a 2 

where Jo is a modified Bessel function of the first kind. The four unknown 
parameters (A, B, C, D) are obtained using the modified HIJING samples. 

iii) Determination of the mean and variance of the true vi distribution, /(f 2 ): 
To extract the fluctuations, we assume a Gaussian distribution for ffa), with two 
parameters, (u 2 ) and a V2 . For given values of (v 2 ) and a V2 , it is possible to take the 
integral in Equation (1) to obtain the expected distribution, g e 7L P (v% hs ; (ify), o~ V2 ). By 
comparing the expected and observed distributions, the values for (i> 2 ) and o~ V 2 are 
found by a maximum likelihood fit. 

The analysis chain is applied to the a/snn = 200 GeV Au+Au data set from 
Run 4 in bins of centrality. The results ((i^) and a V2 ) are obtained separately for 
triangular and trapezoidal V2(n) shape and averaged over bins of collision vertex position. 
The systematic errors (all sources added in quadrature) are estimated by including 



Elliptic flow fluctuations in ^s^n — 200 GeV Au+Au collisions at RHIC 



4 



0.1 



0.08 



cm 0.06 



0.04 



0.02 



PHOBOS preliminary 

Au+Au, 200 GeV, hf, (n|<1 



Event-by-event 
Hit-based 
Track-based 




50 100 



150 200 
Npart 



250 300 350 



Figure 1: Mean flow, (v 2 ), as a function of 
centrality, for ^/snn = 200 GeV collisions at 
mid-rapidity, measured by the event-by-event 
analysis, compared to the published results 
obtained with the event-averaged, sub-event 
based technique using hits or tracks [7] . 
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Figure 2: Relative flow fluctuations, a v% j {v 2 }, 
as a function of centrality, for -^/snn = 200 
GeV collisions at mid-rapidity, compared to 
the prediction, a e J (e), from the participant 
eccentricity, and to an estimate of iV par t- 
induced fluctuations using a fit of (v 2 ) (iVpart)- 



variations from different vertex and (J)q bins and changes introduced by the triangular 
and trapezoidal 1)2(1]) shapes. Since the functional form of the true distribution is 
unknown, differences arising from a flat rather than a Gaussian ansatz for f(v 2 ) are 
included. Furthermore, we have performed extensive studies of the analysis response 
to MC samples, prepared to match the observed multiplicity distribution and {v 2 ), 
in bins of centrality and known input values of fluctuations (including zero). The 
discrepancies of the obtained to the input fluctuations evaluated per centrality bin are 
an additional source of the systematic error. Their contribution becomes increasingly 
important for low values of v 2 , and therefore the most central collision bin (0-6%) 
was excluded from the analysis. Fig. 1 shows the mean flow results obtained with 
the event-by-event analysis projected to mid-rapidity by (v 2 ) = 0.5 (vf 1 ) + ( u 2 rap ))' 
which — within the systematic errors — agree very well with the published results of the 
event-averaged sub-event based technique using hits or tracks [7]. This constitutes an 
important, completely independent verification of the event-by-event analysis. Fig. 2 
presents the relative flow fluctuations, a V2 j (v 2 ). In the ratio, most of the aforementioned 
systematic errors scale out. We find large relative fluctuations of about 40-50%, with 
little centrality dependence. The data are well described by the fluctuations in the 
participants eccentricity, o~ e / (e), obtained in MCG simulations. These are shown as the 
shaded band in Fig. 2, corresponding to 90% confidence level boundaries obtained by 
varying the Glauber parameters within reasonable limits (see Ref. [3]). The contribution 
from iVpart fluctuations is estimated using a fit to the (v 2 ) data folded with the known 
iVpart distributions from the PHOBOS centrality trigger studies. Over the range of 
measured centrality, their contribution (taken in quadrature) is small, less than 15%. 
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3. Conclusions 

First results on event-by-event elliptic flow fluctuations in Au+Au collisions at 
v /5nn = 200 GeV measured with the PHOBOS detector are presented. Large relative 
fluctuations of about 40-50% are found, independent of centrality. The new data are 
well described by fluctuations in the shape of the initial collision region, as predicted 
with the participant eccentricity using MCG simulations. These results substantiate our 
conclusions from previous studies on the relevance of such fluctuations for the elliptic 
flow across nuclear species: The initial- geometry seems to drive the hydrodynamic 
evolution of the system, not only on average, but event-by-event. 
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